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Abstract. In all the microquasars with two hHz QPOs, the ratio of the frequencies is 3:2, supporting our suggestion that a 
non-linear resonance between two modes of oscillation in the accretion disk plays a role in exciting the observed modulations 
of the X-ray flux. We discuss the evidence in favor of this interpretation, and we relate the black hole spin to the frequencies 
expected for various types of resonances that may occur in nearly Keplerian disks in strong gravity. For those microquasars 
where the mass of the central X-ray source is known, the black hole spin can be deduced from a comparison of the observed 
and expected frequencies. 
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INTRODUCTION 

Several Galactic low-mass X-ray binaries (LMXBs) ex- 
hibit quasi-periodic variability (QPOs) of their X-ray 
fluxes, with pairs of « 1kHz frequencies typical in 
the neutron-star sources |JJ] . Kluzniak and Abramowicz 
suggested that these twin kHz QPOs are a mani- 
festation of non-linear resonance that can occur between 
modes of oscillation in an accretion disk in strong field 
Einstein's gravity, but not in Newton's 1 jr potential, and 
pointed out that on this hypothesis the same resonances 
should be present in black hole systems. Pairs of high fre- 
quency QPOs should have been present where only sin- 
gle hHz QPOs had been reported in microquasars. Such 
pairs of several hHz frequencies have now been reported 
in four or five black-hole systems, all in a 3:2 ratio |4], 
substantially strengthening the case for resonance. 

Our originally suggested 01 explanation for kHz 
QPOs in neutron stars was based on these general proper- 
ties of non-linear resonance, which seemed to us to cor- 
respond to the essential features of the observed twin fre- 
quency peaks: 

1. The frequencies of non-linear oscillations v de- 
pend on amplitude, and for this reason they may 
be time dependent and may differ from the fixed 
eigenvalue frequencies v(0) of the system, v(f) = 
v(0) + 5v(r). 

2. Non-linear resonance may occur over a wide fre- 
quency range 5 v. 

3. Both (resonant) frequencies increase or decrease "in 
step" with each other. 

4. The eigenfrequencies of resonant modes are ap- 
proximately in the ratio of small integers, e.g., 2:1. 



These ideas initially received a cool reception (see au- 
thor's note in |3J|), because it was not generally appreci- 
ated that frequency ratios close to 3:2 actually occur 1 for 
kHz QPOs in neutron star sources, and no evidence for 
twin peaks in black hole sources was known at the time. 

In addition to the high frequency QPOs, features in the 
power density spectra can be identified at lower frequen- 
cies, and at least one frequency Vi ow was l° n g known 
to be correlated with V n i g h, one of the kHz/hHz fre- 
quencies in neutron-star/black-hole systems |6]. Quasi- 
periodic modulations of the flux (dwarf nova oscilla- 
tions, DNOs) were first discovered in cataclysmic vari- 
ables (white dwarfs) and these are analogous in many 
respects to the QPOs in LMXBs |7]. However, unlike 
the kHz/hHz QPOs, the highest frequency DNOs do 
not come in pairs. This is consistent with the idea that 
the high frequency pairs arise in accretion disks only in 
strong gravity. 

The relativistic resonance model of black hole QPOs is 
based on fundamental features of strong gravity. Today, 
it is motivated by observations that sharply illuminate the 
physical nature of QPOs: 

1. The correlation V] ow = 0.08Vhi g h between low and 
high frequency QPOs in black hole, neutron stars, 
and white dwarf sources extending over six orders 
of magnitude [H El, proves that in general the 
QPOs are a hydrodynamic phenomenon, and can- 
not be attributed to mere kinematic effects, such 
as Doppler modulation of emission from isolated 
bright spots. (Vi ow may be the "ninth wave" 1 10].) 



Indeed, our paper on this appeared in print only much delayed [5]. 



2. The frequencies of twin peak hHz QPOs in micro- 
quasars seem to scale with mass |4], v ~ 1 /M (Fig- 
ure 1). If true, this would prove their relativistic ori- 
gin. 

3. In all four microquasars with twin peak hHz QPO 
pairs, v U pper/ Viewer = 3/2 (Table 1). 

Although suggestive of a resonance, the ratio 3/2 could 
also be a signature of overtones (flute modes) full , or of 
higher modes of an MHD instability fl2ll at a 'transi- 
tion radius' r* in the innermost part of the disk, which 
excites quasi periodic oscillations with mode frequen- 
cies V ~ n Vk(/*) (in contrast with the resonance model, 
neither of these two models predicts 1 /M scaling, with- 
out making some ad hoc assumptions). However, there 
are additional properties of non-linear resonances which 
may help in their identification. In a non-linear resonance 
combination frequencies, e.g., V U p pe r ± Vi ower , and sub- 
harmonic frequencies may be present 0, e.g., Vi OW er/2- 

Our resonance model may also be applied to twin peak 
QPO sources in neutron stars 1 13]. In refs. 0111, and 
in these Proceedings 1 16], we discuss a resonance in an 
accretion disk or torus excited by an external forcing 
by a millisecond pulsar. A similar forcing is crucially 
important in a different, non-relativistic resonance model 
suggested by Titarchuk 1 17]. 



TABLE 1. The four microquasars in which 
two hHz QPOs are observed. They all have 3:2 
ratio of frequencies. Source of data: |4, 18] 
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FIGURE 1. The IJM scaling of the pairs of QPOs with the 
3:2 frequency ratio |4]. The upper frequency is shown. 



THE ORBITAL AND EPICYCLIC 
MOTIONS. 

Consider a black hole 2 with the mass M$ and angular 
momentum Jo- Inside thin, almost Keplerian accretion 
disks, matter spirals down the central black hole along 
stream lines that are located almost on the equatorial 
plane 6 = 6q = 7r/2, and that locally differ only slightly 
from a family of concentric circles r = ro = const. The 
small deviations, 8r = r — ro, 80 = 9 — 0q are governed, 
with accuracy to linear terms, by 

8r+(0^8r= 8a r , S6 + (Og 1 S0 = Sag. (1) 

Here, the dot denotes a time derivative. For purely Ke- 
plerian (free) motion 8a r — 0, dag — and the above 
equations describe two uncoupled harmonic oscillators 
with the eigenfrequencies (Og = 2nvg, co r = 2nv r shown 
together with the Keplerian orbital frequency, Q. = 2tzVk, 
in Figure 2 for a non-rotating black hole, and in Figure 3 
for a moderately rotating one. 



2 We rescale mass with M = GMq/c 1 = ro, angular momentum with 
a = Jqc/(MqG). We use Boyer-Lindquist coordinates, t,r,6,(j), and 
rescale the radius with x = r/ ra- 



in Newton's theory with the —GMojr potential, 
{2kY 1 GMqi V 3 / 2 = Vk = v r = Ve, but in the strong 
gravity of a rotating black hole, for orbits of the same 
sense of rotation, Vk > Vg > V r . The radial epicyclic 
frequency v r goes to zero at the Innermost Stable Cir- 
cular Orbit for the Keplerian (free) motion, and has a 
maximum at a particular circular orbit outside the ISCO, 
its location depends on the black hole spin | lift 12(1 121I1 . 



1 /M SCALING 

Before the RXTE era, Kluzniak, Michelson, Wagoner 
12211 suggested that the orbital frequency close to the 
marginally stable orbit may be directly observed as a 
QPO, once instruments with sufficiently high time res- 
olution are built, and pointed out that the frequency will 
be inversely proportional to the mass of the compact ob- 
ject. The latter statement applies to any characteristic fre- 
quency in general relativity. 

For example, note that for black holes, all three orbital 
frequencies: Keplerian Vk, radial epicyclic v, , and verti- 
cal epicyclic Vg , also have the general form 




FIGURE 2. Orbital v K , epicyclic radial v r , and vertical Vg 
frequencies for a non-rotating black hole (spin a = 0) with the 
mass Mo = IOMq. Location of the photon radius is marked by 
r p |j, and that of the horizon by rjj. There are no orbits for r < r p j,. 




FIGURE 3. The same as in Figure 2, but for a 1OM black 
hole with a moderately high spin, a = 0.8. 



with a the dimensionless angular momentum of the black 
hole, and f(x,a) a dimensionless function, different for 
each frequency. For all relativistic frequencies, x = x(a) 
is fixed, and then the above formula predicts that frequen- 
cies scale as v = (1 /M)F(a). In particular, each orbital 
resonance n : m discussed later in these contribution oc- 
curs at its own resonance radius x„ :m (a), as shown in 
Figure 13, while at the marginally stable orbit (ISCO) 
V r (r ms ) = 0. All of the models discussed below follow 
the 1 /M scaling. 



NON-RESONANT MODELS 
The highest possible orbital frequency 

In an accretion disk, matter moves, roughly, on circular 
orbits in the region r > r,„ and free-falls in the region 



r < ri„. The transition radius r = ri„, closely coinciding 
with the sonic point, is often called the inner radius of 
the accretion disk. Thin, standard Shakura Sunyaev disks 
with their high efficiencies have an inner edge located 
almost exactly at ISCO. In general, the inner edge is 
located between the ISCO and the marginally bound 
(RISCO, at r m b) circular orbits 12311 . depending on the 
disk efficiency. ADAFs, with their very low efficiencies, 
have the inner edges almost exactly at RISCO 12411 . The 
same is true for super-Eddington slim iz^l and thick 
l23ll disks. For a non-rotating black hole one has r m b = 
AM, VtfOmb) = 4O37(M/AT ) -1 [Hz], and r ms = 6M, 
V^(r ms )-2197(M/M (3 )- 1 [Hz]. 
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FIGURE 4. ISCO orbital frequencies compared with the 
observed upper frequency of the hHz pair of QPOs. 
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FIGURE 5. RISCO orbital frequencies compared with the 
observed upper frequency of the hHz pair of QPOs. 

Let us ignore for the moment that the frequency pairs 
are in a 3:2 ratio. Numerous studies have tried to in- 
fer the properties of neutron stars on the assumption 
that the upper kHz QPO frequency is close to the ISCO 
frequency (e.g., |26j). Identifying RISCO frequencies 
with the upper frequency of the hHz pairs in micro- 
quasars, whose mass is known, would require counter- 



rotating disks (Fig. 5). Comparing the same frequency 
with the one at ISCO leads to sensible but small values 
of the black hole spin (Fig. 4). In most models, the QPO 
frequency corresponds to a characteristic frequency not 
very close to the inner edge of the disk. This leads to sub- 
stantially higher values of the black hole rotation rate. 



The trapped modes 

One of the characteristic properties of the oscillations 
of relativistic disks is the presence of trapped mode 

oscillations l2(il EH S3 Si Sill Ell. The physical 

reason for the trapping is that the radial epicyclic fre- 
quency, v r , is not monotonic but has a maximum value, 
Vtrapp, at a radius r^p larger than the ISCO. For the 
non-rotating black hole hole rtrap = 8M 112 ill . The g- 
mode (inertial-gravity) oscillations 1 30 ] can be charac- 
terized by a restoring force that is typically dominated by 
the net gravitational-centrifugal force. The axisymmetric 
(m = 0) g-modes are centered at r tlap . Non-axisymmetric 
trapped g-modes with the azimuthal wave-number m = 1 
have frequencies J28il . 



V ~V K 0trap)± Vtrapp, and V ~ V K (rtrap). (3) 

In Figure 6 we show the highest frequency connected to 
these oscillations, v upp = VK(r tt ap) + Vtrapp, and compare 
it with observations. 
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FIGURE 6. Maximum frequency of the trapped, m = 1, 
g— mode compared with the upper QPO frequency. 



Dragging of inertial frames and the c-mode 



The 'corrugation' c-mode 1311 1271 13211 is a non- 
axisymmetric, vertically incompressible wave near the 
inner edge of the disk that exists only for co-rotating 
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FIGURE 7. Frequency of the c-mode at ISCO compared 
with the upper QPO frequency. 



disks, a > 0. It precesses around the angular momen- 
tum of the black hole. Its frequency coincides with the 
Lense-Thirring frequency produced by the dragging of 
inertial frames. In Figure 7 we compare with observa- 
tions the highest frequency connected to the c-mode, as- 
suming that the mode locates at ISCO: 



K 11PP 



VLT(/ n 



ac r G 3 
7Cr G r ms 



(4) 



In reality, the mode is trapped further out in the disk, and 
correspondingly, the inferred value of black hole spin is 
higher |32]. 



NON-LINEAR, RELATIVISTIC 
ORBITAL RESONANCES 

"Keplerian" resonances 

It is possible for the radial epicyclic frequency to be in a 
resonant relation with the orbital frequency, nv r = mVK, 
with n, m integer 13311 . For example, g-modes have pat- 
tern frequency 2nv,„ = ±(v r ±toVa-), and these can be in 
co-rotation resonance, i.e., with v,„ = Vk l34ll. The case 
Vjf/v r = 3/2 is excluded by observations 1 33 ] (Fig. 8). 
The remaining possibilities 1 33 ] are that the upper fre- 
quency is Vk ± v r , with Vk/Vv = 2, or V/f/v,- = 3, and 
Vupper/Vr = 3/2 (Figs. 9, 10). However, co-rotation res- 
onance leads to damping, and not excitation, of modes 

Another possibility is based on the following idea 
(3(1. When the potential vorticity is conserved, coherent 
vortices tend to form in pairs with opposite vortices. One 
may imagine that because the spatial distance between 



the two structures, which oscillates with the epicyclic 
radial frequency, depends on the velocity profile of the 
disk, i.e., also on the oscillations of orbital velocity, a 
resonance between these two frequencies is possible. 
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FIGURE 8. Unsuccessful fit of observations to the Keplerian 
3:2 resonance. 



Epicyclic motions 

The effective potential U(r,9;£) for orbital motion of 
a particle with a fixed angular momentum I > l ms has 
a minimum at tq{€), corresponding to the location of a 
stable circular orbit. The second order term in its Talyor 
expansion (for simplicity we write it on the equatorial 
plane 9 — 7t/2) gives the epicyclic frequencies, terms in 
the next order 



d 2 U 
dr 2 



dr 3 



( r -r ) 3 + ... (5) 



contain higher than quadratic terms, which means that 
small oscillations around the minimum at r — rr, are 
described by non-linear differential equations 1371 13811 . 
Non-linear resonances that may be excited in these non- 
linear oscillations have several characteristic properties 
that closely resemble those observed in QPOs. 
We consider two possibilities. 
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FIGURE 9. Keplerian 3:1 resonance. 




The forced 3:1 and 2:1 resonances 

A direct resonant forcing of vertical oscillations by 
the radial ones through a pressure coupling, and with 
8r ~ cos(o r f ), is evident in recent numerical simulations 
of oscillations of a perfect fluid torus 11511 . This supports 
a possible model for the twin peak kHz QPOs: a forced 
non-linear oscillator, 



8d + (O g 1 86 + O 2 (80) = hcos((Q r t), (o e ~n(o r . (6) 

Obviously, there is no integer value of n such that (Og and 
(O r could be in the 3:2 ratio. However, non-linear terms 
alow the presence of combination frequencies 1371 13811 . 



00- 



(Og - 0) r , (0 + = (Og + (O r - 



(7) 



One of these combination frequencies may be in a 3:2 
ratio with (Og if and only if n = 2, or n = 3 in this forced 
resonance. Simple arithmetic shows that in these two 
cases the observed frequencies V] ower = (O\ ov/er /2n and 



Vupper — 



Cupper /2 it are uniquely given by, 



[3 : 1] 0*0 



(0 = 2 (0 r , Copper = (0g=3(0 r , (8) 



log(A//M 



[2 : 1] ©upper = C0+ = 3 (O r , (Blower = 00g = 2(O r . (9) 



FIGURE 10. Keplerian 2:1 resonance. 



We fit observed QPOs to these predicted by the forced 
epicyclic 3:1 and 2:1 resonances in Figures 11 and 12. 



M/M Q 
10 



12 14 ie 




0.9 1 1.1 

log (M/M a ) 







•1 c 
1 D 






14 


oUUnz 




- 450Hz 




I Z 


- 400Hz 




- 350Hz 






- 300Hz 




10 








- 250Hz 










8 


- 200Hz 










6 


- 0.608 










4 


- 0.457 






- 0.305 a 






- 0.154 




2 


- 0.002 






- -0.150 












1.3 



FIGURE 11. Fit of the predictions of the 3:1 epicyclic forced 
resonance model to observations. 
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FIGURE 12. The 2: 1 forced epicyclic resonance. 



3:2 parametric resonance 

We shall start with an argument appealing to physical 
intuition and showing that the resonance to be discussed 
now is a very natural, indeed necessary, consequence of 
strong gravity. In thin disks, random fluctuations have 
8r 3> 89. Thus, 8r80 is a first order term in 86 and 
should be included in the first order equation for vertical 
oscillations Q. The equation now takes the form, 



S8 + a>l[\+h8r]Se = 8a e , 



(10) 



where h is a known constant. The first order equation for 
8r has the solution 5r = Aocos(o r f ). Inserting this in 
( ITOl together with 8a g = 0, one arrives at the Mathieu 
equation (Aq is absorbed in h), 



89 



colli 



-hcos(co r t)}89 = 0, (11) 

that describes the parametric resonance. From the theory 
of the Mathieu equation one knows that when 
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FIGURE 13. The radii of three epicyclic resonances, 3:2, 3:1 
and 2:1. Also shown are the location of the centre of largest 
possible accretion torus (dashed lines), and the location (dotted 
lines), closer to the horizon, of the maximal locally emitted flux 
from the standard thin disk. 
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(12) 



the parametric resonance is excited l38ll . The res- 
onance is strongest for the smallest possible value of 
n. Because near black holes v, < Vg, the smallest pos- 
sible value for resonance is n = 3, which means that 
2 Vg = 3 v r . This explains [39] the observed 3:2 ratio, if, 
obviously, 



Vupper — Ve, Vi 0W er — V r . 



(13) 



Parametric resonance of the type discussed above was 
found in numerical simulations of oscillations in a nearly 
Keplerian accretion disk by Abramowicz et al. 14011 and 
confirmed by exact analytic solutions num. The ana- 
lytic solution is accurate up to third order terms in 8r, 
86, and based on the method of multiple scales 1 38] . 
Existence of the 3:2 parametric resonance is therefore a 
mathematical poroperty of thin, nearly Keplerian disks. 

It was found that the resonance is exited only in the 
non-Keplerian case, with some weak forces Sag ^ and 
8a r ^ present. Their origin is certainly connected to 
stresses (pressure, magnetic field, viscosity), but how ex- 
act details remain to be determined — at present Sag 
and 8a r are not calculated from first principles but de- 
scribed by an ansatz 3 . Of course in real disks neither 



While the lack of a full physical understanding is obviously not 
satisfactory, the experience tells that such a situation is not uncom- 
mon for non-linear systems. Examples are known of mathematically 
possible resonances causing damage in bridges, areoplane wings etc., 
for which no specific physical excitation mechanism could have been 
pinned down 1331 . 
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FIGURE 14. Fit of the frequency V upp = Vg (M.a) predicted 
by the 3:2 parametric resonance model to the frequencies ob- 
served in three microquasars with known masses. The spin pa- 
rameter a is not known from observations, and the lines a = 
const are calculated from the model. One should notice that the 
deduced black hole spins are rather high, and fall in a rather 
narrow range. 
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FIGURE 15. The 1/M scaling of the frequencies of 3 : 2 QPO 
pairs may be used for mass estimates (from Abramowicz, Kluz- 
niak, McClintock, & Remilard, 2003). 



8r = Aocos(ovf), nor Sag = exactly, but one may 
expect that because these equations are approximately 
obeyed for thin disks, the parametric resonance will also 
be excited in realistic situations. And this is indeed the 
case 0. 

The parametric resonance occurs at a particular ra- 
dius 7"3 : 2(a), determined by the condition 3co r (ri-2,a) = 
2a)e(r3 2,fl). We show the function r^ia) in Figure 13. 
In Figure [H] we fit the 3:2 resonance theoretically pre- 
dicted frequencies to the observational data for the three 
microquasars with known masses. The scatter for the par- 
ticular 3:2 resonance is not very large because (Fig. 13) 
this resonances occurs at xj-2{a) > 4, where the influence 
of a is not dominant. 



TABLE 2. Black hole spin in three micro- 
quasars calculated by fitting observations to the 
three resonance models 



Microquasar 


3:2 


2:1 


3:1 


XTE 1550-564 


0.94 


0.27 


0.46 


GRO 1655-40 


0.96 


0.36 


0.55 


GRS 1915+105 


0.84 


0.02 


0.23 



APPLICATIONS 

The 1/M scaling of the twin peak QPOs frequencies with 
the 3:2 ratio, was proposed by Abramowicz, Kluzniak, 
McClintock & Remillard |44] as a method for estimat- 
ing black hole masses in AGNs and the recently dis- 
covered ultraluminous X-ray sources (ULXs), based on 



Mirabel's analogy between microquasars in our Galaxy 
and distant quasars 1 4-5 ] . Indeed, if the analogy is also 
valid for accretion disk oscillations, then discovering in 
ULXs the twin peak QPOs frequencies with the 3:2 ratio, 
would resolve the controversy about their mass: if ULXs 
black holes have the same masses as microquasars, the 
frequencies will be ~ 100 Hz; if ULXs black holes are 
~ 1000 times more massive, the frequencies will be 
~ 0.1 Hz instead. 
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